Slowing down Josephson vortex lattice in Bi2Sr2CaCu208+5 with pancake vortices 
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We study theoretically and experimentally influence of pancake vortices on motion of the Joseph- 
son vortex lattice in layered high-temperature superconductors. Mobility of the Josephson vortices 
in layered superconductors is strongly suppressed by small amount of pancake- vortex stacks. Mov- 
ing Josephson vortex lattice forces oscillating zigzag deformation of the pancake- vortex stacks con- 
tributing to damping. The salient feature of this contribution is its nonmonotonic dependence on 
the lattice velocity and the corresponding voltage. Maximum pancake effect is realized when the 
Josephson frequency matches the relaxation frequency of the stacks. The pancake- vortex damping 
is strongly suppressed by thermal fluctuations of the pancake vortices. This theoretical picture was 
qualitatively confirmed by experiments on two mesas prepared out of Bi2Sr2CaCu2 08+5 whiskers. 
We found that the Josephson- vortex flux-flow voltage is very sensitive to small c-axis magnetic field. 
The pancake-vortex contribution to the current indeed nonmonotonically depends on voltage and 
decreases with increasing temperature and in-plane magnetic field. We also found that irradiation 
with heavy ions has no noticeable direct influence on motion of the Josephson vortices but dramati- 
cally reduces the pancake- vortex contribution to the damping of the Josephson vortex lattice at low 
temperatures. 



I. INTRODUCTION 

The layered crystalline structure of the cuprate high- 
temperature superconductors leads to existence of two 
types of vortices in these materials, pancake-vortex (PV) 
stacks^ induced by the magnetic field component per- 
pendicular to the layers and Josephson vortices (JVs)^ 
created by the magnetic field component parallel to the 
layers. Repulsive interaction between the vortices of each 
type results in formation of regular vortex lattices. In 
particular, the magnetic field applied along the layers 
generates triangular lattice of the JVs stretched along 
the layer direction. The anisotropy factor 7 sets the im- 
portant field scale, Bcr = <^>o/(27r7s^) , where s is the 
interlayer periodicity. When the magnetic field exceeds 
Bcr the Josephson vortices homogeneously fill all layers 
forming a dense lattice^. In highly anisotropic materials, 
like Bi2Sr2CaCu208+5 (BSCCO) this field scale is rather 
moderate ~ 0.5 tesla. 

In BSCCO, due to a very weak Josephson interlayer 
coupling, two types of vortices can coexist in the tilted 
magnetic field^. The static attractive interaction be- 
tween JVs and PV stacks^ leads to many exotic vortex 
states, such as mixed chain-lattice state^&i and pancake- 
chains state^, and has been used for visualization of JV 



stacks 
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Dynamic properties of the crossing-lattices state have 
been studied in much less details. A particularly inter- 
esting topic is dynamic behavior of the JV lattice. An 
external transport current flowing across the layers drives 
the JV lattice along the layers. Due to relatively weak 



intrinsic dissipation, the Josephson vortices can be ac- 
celerated up to very high velocities. Dynamics of the 
JV lattice in BSCCO have been extensively studied by 
several experimental groups (see, e.g., Refs. Il2lll3lll4l) . 
When magnetic field is tilted at small angle with respect 
to the layers, the c-axis field component generates small 
concentration of PV stacks. Alternating in-plane super- 
currents of static JV lattice force zigzag deformations of 
the PV stacks^^, see Fig.^ It is well known that mobility 
of JVs is strongly suppressed by a very small amount of 
PV stacks^^'^'''^^. As a consequence, studying the JV lat- 
tice transport always requires a very accurate alignment 
of the magnetic field with the layers. In spite of that 
common knowledge, JV lattice dynamics in presence of 
the PV stacks has never been investigated systematically. 

In the case of strong anisotropy, the JV lattice can 
move through static PV stacks. Even in this case the PV 
stacks will strongly impede motion of the JV lattice. Dy- 
namic behavior of the PV stack naturally follows from its 
static configuration. The moving JV lattice forces oscil- 
lations of the PV stacks leading to additional dissipation 
and slowing down the lattice. In this paper we investigate 
this effect quantitatively in the case of dense JV lattice. 
Influence of the PV stacks on motion of an isolated JV 
has been considered theoretically in Ref. iIQl 

The paper is organized as follows. In Sec. ^ we present 
theoretical models describing influence of the PV stacks 
on motion of the dense JV lattice. We compute the dy- 
namic friction force generated by PV stacks and study 
suppression of this force by PV fluctuations. We also 
consider influence of the PV fluctuations on the regular 
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FIG. 1: Sketch of a pancake- vortex stack in presence of the 
dense JV lattice. The arrows illustrate current distribution 
of the JV lattice and the discs illustrate pancake vortices. 
The JV lattice induces zigzag deformation of the stack. Mov- 
ing lattice induces oscillating zigzag deformations which con- 
tribute to dissipation, see animation . 



JV flux-flow resistivity and influence of columnar defects 
on the PV-induced damping of the JV lattice. In Sec. 
mil we present experimental results. Studying the flux- 
flow of the JV lattice for small c-axis magnetic fields, we 
extracted the PV contribution to the JV damping and 
studied its dependence on the voltage, temperature, and 
in-plane field. We also found that this PV contribution 
is strongly suppressed by heavy-ion irradiation. In Sec. 
IIVI we discuss comparison between the experiment and 
theory and possible applications of the studied effect. 



II. THEORETICAL MODELS 

A. Basic model 

A general description of JV dynamics in presence of the 
PV stacks is rather complicated. We consider first the 
simplest situation, in which influence of the pancake vor- 
tices on motion of the JV lattice can be traced in a clear- 
est way. We assume that a strong magnetic field, B > 1 
tesla, is applied in the x-z plane at a very small angle with 
respect to the layers. Such a magnetic field generates the 
dense JV lattice and a dilute array of the PV stacks. We 
neglect both pinning and thermal fluctuations of the pan- 
cake vortices. We start with derivation of the interaction 
between the JV lattice and deformed PV stack via the 
inter layer Josephson energy, £j = —Ej^^J (Pr cos On 
where On is the phase difference between the layers n 
and n -\- 1. The static dense JV lattice is characterized 
by the phase difference, Ojn{y) ~ kniv — yj) + ttu with 
kn = 27TsBx/^o and yj describing the JV lattice dis- 
placement. The deformed PV stack with the coordinates 

Un {U 

n,x-)'^n,y) generates the additional phase differ- 
ence 



6>p(r, Un+i, Un) = arctan 



y-Un+l,y 
X '^n+l,x 



arctan 



y-ur^ 



and modifies the Josephson energy. In addition to inter- 
actions mediated by interlayer Josephson coupling, PVs 
are also aligned by the magnetic interaction^. Due to its 
long-range nature, this interaction can be well approxi- 
mated by a simple parabolic potential (see, e.g., Ref. ^3) 
with the magnetic-spring constant, Kp ^ s^'^C/ (47rA^) 

where C ~ \n.{\/rw) and = (('^n+i — '^n)^)"'^^^ is the 
typical wandering distance. Therefore in a homogeneous 
superconductor the total energy change per one layer 
due to PV stack deformation is given by the sum of the 
Josephson and magnetic energies 



N ^ 



K 

(-l)"" £;jC(Un+l, Un, yj) + 



, (1) 



where N is the total number of layers and the func- 
tion C(u2, ui, yj) in the lowest order with respect to the 
Josephson coupling energy, Ej^ is given by 

C{u2, ui, yj) = J d^r{cos [kniy-yj)] 

- cos [kH{y - yj) + 6>p(r, U2, ui)]} (2) 



/(/C//U12) 
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with U12 = U2 — ui and 



/(v) = J dxdy ^1 — — 



Vxy - VyX 



|r- v/2| |r + v/2 



-v/2||r + v/2| 
sin 7/ 



cos y 



(3) 



with V = kHy^i2- In the regime kH\y^i2 \ ^ 1, this dimen- 
sionless function can be computed analytically up to the 
third order with respect to the reduced variable v. 



7(v). 
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Vy and Ch = 8exp(-7E) 
is 



(4) 
^ 4.492 
The 
linear- 



with — 

(7e ~ 0.5772 is the Euler constant), 
linear in Vx term in /(v) gives the 
displacement contribution to the Josephson energy, 
- {47rEj/NkH)^ni~'^T ^n,xCos[kHyj]' This term 
describes forces acting on the straight PV stack from 
the alternating in-plane currents induced by the JV 
lattice, jy{y,n) ^ (-1)'^ jh cos [kniy - yj)], where 
jh = (27//i)jj, h = 27r5a;75^/$o is the reduced 
magnetic field and jj = {27tc/^o)Ej is the Josephson 
current density. Due to this term the ground state 
corresponds to the alternating PV deformations along 
the direction of the in-plane field^^, Un,x = (— l)^'^a 
with Ua = {s^o/c)jh/Kp = 4X'^/[hjs\n{X/ua)], see Fig. 
n For such deformations, the quadratic term in the 
Josephson energy cancels out and does not influence 
the deformation amplitude. The assumed condition 
knUa < 1 is satisfied if 7 > A/5. 
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Lets consider now dynamic behavior. A transport 
current applied across the layers drives the JV lattice 
along the layers. We consider the lattice slowly mov- 
ing with constant velocity through the PV stacks along 
the y axis, yj = vt. Such motion generates the elec- 
tric field, Ez = Bxv/c. The lattice velocity is assumed 
to be much smaller than the Swihart velocity so that 
the lattice preserves its static structure. The homoge- 
neously moving JV lattice forces PVs to oscillate with 
the Josephson frequency, Uxiji^t) = (— l)^iXa(t), with 
Ua{t) = Uao<^os{iOEt + a) and coe = knv. For a PV 
stack located near y = 0^ the alternating deformation, 
Uaif)^ is described by a simple oscillator equation 



r]pUa + KpUa = fh cos , 



(5) 



where fh = {s^o/c)jh = Air^sEj/h and r]p is the pan- 
cake viscosity coefficient. Solution of this equation is 
given by 



Uait) Re 



fhex.p{iujEt) 



iujeVp 



(6) 



The frequency response of the stacks is determined by 
the relaxation frequency, co^ = Kp/rjp. The oscillation 
amplitude drops at cc;^; > ujr- 

The average friction force per PV acting on the moving 
JV lattice is given by 



d5£ \ 
dyj / 



Nk 



^— ^ ({-'^r ^ [^H (Un+l-Un)] siu 



UOEt-kH 



(7) 



where (...) notates time averaging. In the regime when 
only small alternating deformations in the x direction are 
present, this gives J^y = — 47r£^j (^/^(t) sin (cc;^;^)), and, 
using the solution Q, we obtain 



J^y — —27rEjfh 



^eVp 
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The average velocity of the JV lattice is connected with 
the applied current, j^, by the force-balance condition 



{Bjc)jz = 



(9) 



where r]jf[ is the bare flux-flow friction coefficient of the 
JV lattice and Up = B^/s^q is the concentration of PVs, 
and the third term is the driving force from the trans- 
port current. This corresponds to the following current- 
voltage dependence 



3z{Ez)=cJmEz^5j{Ez) 



(10) 



where crjff = (cjBx) T^jff is the bare flux- flow conductiv- 
ity and 



(11) 



is the current enhancement due to the PV-induced damp- 
ing of JV motion, which we will call excess pancake- vortex 
current^ with cjp = ap{Bz^Bx) = 27r'^npr]ps'^j'^/Kp and 
Vr = [^o/ {27Tc)]Kp/r]p is the voltage drop per junction 
corresponding to the relaxation frequency. The electric- 
field dependence of the conductivity Sj{Ez)/Ez resem- 
bles a well-known Drude frequency dependence. In- 
troducing the pancake flux- flow conductivity a{^{Bz) = 



ripC^ / {s^qBz)^ we can express the conductivity and volt- 
age scales, (jp and F^, via experimentally-accessible pa- 
rameters 



2 ' 



Vr 



C*0 



87r2A2<Tff(BA) 



with Bx = <l>o£/(47rA^). Note that dp scales with the 
field components as ap oc Bz/B^. 

The key feature of the excess PV current, dj{Ez)^ is 
that it depends nonmonotonically on the electric field. 
The maximum current due to the pancake vortices can 
be estimated as jmax = jj{Bz/B\) jln?. As Gp oc Bz, 
the total I-V dependence ([Tn|) becomes nonmonotonic at 
sufficiently large Bz. In fixed-current experiments this 
leads to voltage jumps. Eq. ([Tn|) also determines the 
angular dependence of voltage at fixed current frequently 
measured experimentally. Introducing the tilt angle of 
field, ^ Bz/Bx <C 1, we can rewrite Eqs. (tTHj) and (tTT|) 
in the form 



= h(l^ {sEz/Vrf 



jz - crmEz 
aiEz 



(12) 



with (Ji = 4<jff (5a)A^/[7^5^£]. This equation explicitly 
determines the dependence Ez{0) at fixed jz. In par- 
ticular, in the case Vr <C sj^/crjff, the maximum and 
minimum of the dependence 0{Ez) are given by 

^^ 

at Ez ^ jz/'^crm 



h {sjz 



2hsj 



aiVr 



at Ez ~ Vr/s 



At these angles the voltage will have jumps. As we see, 
in a simple model both angles are proportional to the in- 
plane field. The angular dependence of the electric field 
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FIG. 2: The angular dependence of the electric field at fixed 
current following from Eq. (fT^ at Vr <^ sjz/ajf^. Voltage 
jumps are expected at the angles ^max and ^min 



following from Eq. (|12|) at Vr <C sjz /(Jjf[ is illustrated in 
Fig.H 



B. Friction force from fluctuating pancake- vortex 
stacks 



In real experimental conditions the result may be 
strongly influenced by pinning and thermal fluctuations. 
Let us consider influence of thermal fluctuations. At a fi- 
nite temperature, in addition to the regular zigzag defor- 
mations, the PV stack has random fluctuating displace- 
ment Uf^n = {uf^x,n^Uf,y,n)' Thcsc paucakc fluctuations 
influence the damping of the JV lattice in several dif- 
ferent ways. The most direct channel is that, due to 
random displacement in the y direction (perpendicular 
to the field), exact time periodicity of the force from 
the JV lattice acting on PV is perturbed. This leads 
to the Debye- Waller- like suppression of the pancake os- 
cillation amplitude and reduction of the corresponding 
current enhancement. The PV displacements also mod- 
ify the interaction force between PVs and JV lattice (in 
the low orders with respect to displacements this effect is 
described by Eq. Q)). Both above effects have the same 
order and become strong when the typical PV fluctuation 
becomes of the order of the JV lattice period. The pan- 
cake fluctuations also perturb phase distribution of the 
JV lattice leading to renormalization of its bare flux-flow 
resistivity. We will also estimate this effect. 

We now analyze the fluctuation effects quantitatively. 
At finite temperature and for the homogeneously mov- 
ing JV lattice, yj = vt, the PV displacements obey the 
following equation 



r]pUn + KpUn- 



(13) 



where the function C(ui,U2,^j) is defined by Eqs. (O 
and and fL,n(^) is the thermal Langevin force with 



the correlation function 

This equation generalizes Eq. Q for finite temperatures 
and large pancake displacements. Introducing the di- 
mensionless variables for the time, t = cort and coordi- 
nate u = knu^ we conclude that the overall behavior is 
determined by the two dimensionless parameters. 



aj 



ttEj 



(14) 



giving the relative strength of magnetic and Josephson 
couplings between PVs in different layers and the reduced 
temperature 



T = 



{2ttsXB, 



f (47rA)' k^T 



(15) 



This reduced temperature has a transparent physical 
meaning, at T ~ 1 the fluctuating pancake displacement 
becomes of the order of the JV lattice period. We can see 
that the relative strength of fluctuations rapidly increases 
with increasing in-plane field. 

We numerically solved the dynamic Langevin equa- 
tions for the pancake displacements (fT3j) with numer- 
ically computed dependence /(v) from Eq. and used 
obtained solution to compute the average friction force 
Q. Figure [3^ shows the Josephson- frequency depen- 
dences of the friction force for different reduced tem- 
peratures for typical value of the parameter aj, aj = 0.1. 
We can see that the main effect of fluctuations is sup- 
pression of the maximum of the J-y{ooE) dependence cor- 
responding to the maximum of 5j{Ez). To quantify this 
suppression, we show in Fig.|3l3 temperature dependences 
of the maximum friction force for two values of the pa- 
rameter a J, 0.1 and 0.2. The maximum force decreases 
by factor of 2 at T ^ 0.5 — 0.7. The maximum only 
slightly shifts to higher frequencies with increasing tem- 
perature. Also, at high temperatures the force decays 
slower than I/loe at large frequencies. 



C. Reduction of the regular flux-flow conductivity 
by pancake- vortex fluctuations 

Pancake fluctuations also perturb regular phase distri- 
bution in the moving JV lattice leading to modification 
of the regular flux-flow conductivity. This gives an addi- 
tional channel of P V influence on the JV mobility. In this 
section we estimate this effect quantitatively. Dynamic 
behavior of the JV lattice in layered superconductors is 
described by the coupled reduced equations for the inter- 
layer phase differences On and reduced fields hn, see, e.g., 
Ref.ill 

^ ' - ' ' - (16a) 







dOn 

dy 



dx dy 
d I dOn hx^n 



dt V dy /2 
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FIG. 3: (a) Evolution of the friction force vs Josephson fre- 
quency dependence with increasing reduced temperature (fT^ 
for aj = 0.1. In this plot Ty is the friction force induced 
by one fluctuating PV stack on JV lattice defined by Eq. Q 
in units of Kp/kn- The frequency unit is Ur. The curves 
are calculated using pancake displacements obtained by di- 
rect numerical simulation of the Langevin dynamic equations 
(|13|) . The upper solid line represents the analytical result Q 
valid at T = and aj <^ 1. (b)The temperature dependence 
of the maximum friction force for two values of the parameter 
aj, 0.1 and 0.2. 



In these equations the units of space and time are given 
by the Josephson length, Aj, and the inverse plasma fre- 
quency, l/cjp, Uc and Vab are the dissipation parameters 
related to the quasiparticle conductivities, and / = Xab/s. 
For the moving triangular lattice the phase differences 
can be represented as 

6>n(r, t) = ujEt - kny + Trn + 6>^,n(r, t) + 0n(r, t) 

where i9^,n(r,t) = 6>p(r - R_^-, Un+i,j (t), Un,j (t)) is 
the phase perturbation induced by PV displacements, 
Un,j(t), of the stack located at Rj , (/)n(r, t) is the correc- 
tion due to the Josephson coupling. Solving Eqs. ([TB|) by 
Fourier transform in the first order with respect to sin^^, 
we obtain 

[sm{u;Et-kHy^7rn^0y^n{r,t))]^^ 
0(k,^,cj)= -— ^ (17) 



with A(k, g', cj) =uj'^ — ivc^ 



A(k,g,cj) 

2(l-COS^) + (l+Zi/a6^)//2' 



where [• iiotates Fourier transformation with re- 

spect to coordinate and time. To obtain I-V character- 
istic we have to evaluate the average reduced Josephson 
current ij = (sin On)^ which in the first order with respect 
to (j)n{Y-,t) is given by 

ij ^ (cos {iJEt - kny + Trn + 6>^,n(r, t)) ^n(r, t)) . 

Note that the PV excess current derived in the previous 
sections also can be obtained using this approach if we 
add to t) the dynamic phase perturbation due to 

the oscillating PVs. 

Using result (|T7)) . we derive 



ij = S{k,q,ij)ijo{kH-K7T-q,ujE-uj) 

J k,g,a; 



(18) 



where 5'(k, g,cc;) is the Fourier transform of the phase 
correlation function 



Sn-n'{r - v',t) = (exp [i {Oy^nir^t) - ^^,n'(r^O))]) 

^(^,,n(r,t)-^,,nKr ,0))')/2 



and 



exp 



ijo(k,g,cj) = 2 



1 



A(k,^,cj) 



In particular, without PV fluctuations 6'(k, g,cc;) = 
S{]i.)S{q)S{uj) and the average Josephson current is sim- 
ply given by ijo(kH^7r^u;E) giving at small coe the bare 
JV flux- flow resistivity (see Ref. 21). 

Let us evaluate the phase correlation function 6'n(r, t). 
Neglecting correlations between displacements in differ- 
ent stacks, we obtain 

Sn{r-r\t) = 1-^ /dR(l-exp [-G(r-R, r'-R)]) 

(19) 

G(r, rO = I ([^p(r, Un+i(t), Un{t))-Op{r\ ui(0), uo(0))]' 
Using expansion of the PV phase difference at r |wn| 

with Wn Un+l-Un, 6>p (r, Un+1 , Un) ^ - [r X Wn]^/r^, 

we can see that the integral over R logarithmically di- 
verges at large R. Cutting off this divergency R = Rc^ 
we obtain with logarithmical accuracy the following re- 
sult 

Sn (r, t) = l- ^n, f{rl)ln^- (w, (t) wq (0)) In ^ 

(20) 

with (r^) = (w^). From Eqs. (|TKj) and we can con- 
clude that the PV fluctuations reduce the regular contri- 
bution to the Josephson current and the relative reduc- 
tion can be estimated as Sjr/j ^ — Trn^ (r^) In (i^c/^w)- 
This correction roughly corresponds to suppression of 
the Josephson coupling in the area 7rr^ around the PV 
stack. At small velocities this reduction is typically much 
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smaller than the current increase due to pancake oscil- 
lations. However, as the PV excess current decays with 
increasing velocity, the regular contribution may become 
dominating at high velocities. The relative contribution 
of the regular term also increases with increasing tem- 
perature. 



D. Role of columnar defects 

Columnar defects produced by heavy-ion irradiation 
are known to be the most efficient pinners of PVs, see Ref. 
[2^ Moreover, they effectively suppress fluctuation and 
align the PV stacks^^. Experimentally, this alignment 
leads to significant increase of the Josephson Plasma Res- 
onance frequency^^. Therefore we can expect that the 
columnar defects will strongly suppress the PV-induced 
damping of J Vs. Naively, one may think that effect of 
the columnar defects can be described by a simple en- 
hancement of the effecting spring constant in the confin- 
ing parabolic potential. However, a more accurate anal- 
ysis below shows that this is incorrect description. A 
columnar defect produces a very strong local potential 
for PVs, which is poorly described by the parabolic ap- 
proximation. Due to a discrete nature of the PV stack, its 
statistical distribution around the columnar defect con- 
sists of two parts^^: a sharp peak centered at the defect 
corresponding to PV located inside the column and a 
very wide envelope function corresponding to PV located 
outside the column. In such situation, we can expect 
that PVs confined inside the column do not contribute 
to JV damping while PVs outside the column give almost 
the same contribution as free PVs. This means that the 
columnar defect simply reduce the excess current by the 
probability factor to find PV outside the columnar defect. 
Pout, as Sj Pout^j' To estimate this probability, we 
consider a fluctuating pancake vortex near the columnar 
defect, which we model as an insulating disk with radius 
b. The ratio of probability to find the PV outside the col- 
umn. Pout, to probability to find PV inside the column. 
Pin, can be evaluated as 



So 



■ exp 



ev{r) 



(21) 



where ev{r) is the PV energy, measured with respect to 
the ground state corresponding to the PV located in the 
column center. 



ey{r) = Up^ seo In (l - b'^/r 



■Kpr^/2. 



with £q = ($o/47rA)^. The first term. Up ^ sso\ri{b/^) is 
the pinning energy for 6 > ^, the second term is the inter- 
action energy with columnar defect, and the third term is 
the magnetic coupling energy, so is elemental area which 
can be evaluated by analyzing the Gaussian fluctuation 
of the order parameter in the vicinity of vortex core^^. 



So 



\n{Hi) sso ' 



in which ^ is the coherence length. Hi = A/^ is the 
Ginzburg-Landau parameter, and ~ 1. In contrast to 
consideration of Ref. 23, in the case of large in-plane field 
we can neglect Josephson interaction in between PVs in 
neighboring layers and take into account only magnetic 
coupling. The value of the integral is determined 
by the competition between the two factors: the large 
energy cost seo^i^{b/^) for putting PV outside the col- 
umn and large area of integration limited by the shallow 
magnetic-coupling parabolic potential. For estimate, we 
can neglect the long-distance tail 5^0 In (l — 6^ /r^) in the 
pinning potential. In this case we obtain the following es- 
timate 



UA 27rT 



exp ( - — 



soKp 



'' exp 



T ) Cs 



Therefore, the probability to find PV outside the column 
can be estimated as 



1 



1 + exp [seoHb/^)/T]Cs/K.^ 



(22) 



Ratio Pout/^in becomes of the order of 1 at the typical 
depinning temperature, T*, which can be estimated as 



seo ln(6/0 
21n(/^) 



(23) 



The columnar defects strongly reduce the damping due 
to PVs at T <C T* and this effect rapidly decreases 
at T > T*. The value of this important tempera- 
ture scale is most sensitive to the London penetration 
depth. For optimally doped BSCCO, assuming A(T) = 
200nm/vTl - (T/Te)^), ^ = 15nm/y(l - (T/Te)^) with 
Tc 90K, and b 7nm, we estimate T* ^ 58K. For 
weakly underdoped BSCCO, taking A(0) = 250nm and 
Tc = 85K, we obtain T* ^ 45K. 



III. EXPERIMENTAL RESULTS 

Josephson flux-flow (JFF) measurements have been 
done on the two mesa structures fabricated from 
the BSCCO single-crystal whiskers using double- 
sided processing by the focused ion beam (FIB) 
technique^^, see Fig. 0] The mesas had geometrical 
sizes LttXL5xLc=26/imx5.8/imx (0.18-0. 23)/im (mesa 
#1) and 20/imx5/imx (0.5-0. 6)/im (mesa #2). The mesa 
#1 was prepared out of the overdoped whisker grown in 
oxygen flow. The mesa #2 was prepared out of the un- 
derdoped whisker grown in a mixture of 10-20 of O2 and 
80-90% of Ar. 

At fixed temperatures and fixed in-plane component of 
the magnetic field {H \\ b) we measured the JV flux- flow 
resistance at several fixed currents and the I-V char- 
acteristics as a function of the c-axis component of the 
field {H II c). The cryostats with two perpendicular coils 
have been used for these experiments, the main coil pro- 
vided the parallel field component, H^^ while another coil 
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FIG. 4: Two mesa structures fabricated out of BSCCO 
whiskers for transport measurements. 



induced the perpendicular component, Hz. The strictly 
parallel orientation has been found using sharp maximum 
of i^ff {Hz ) dependence by sweeping Hz . We estimated an 
accuracy of field rotation in this way to be within 0.01°. 
The experiments have been carried out at high temper- 
atures above 40K to avoid flux trapping effects. As a 
separate part of the experiments, we studied the influ- 
ence of columnar defects on the JFF with and without 
the perpendicular field component. The irradiation was 
made by Pb-ions with energy of 1 GeV at the GANIL 
accelerator (Caen, France). The defects were introduced 
into the mesa #2 along the c axis with concentration 
3 • lO^cm"^ corresponding to the matching field of about 
60 Oe. For a comparison, the JFF characteristics have 
been measured in the mesa #2 first just before irradia- 
tion and then remeasured exactly at the same conditions 
within a week after irradiation. 

Figure [SJleft panel) shows the dependence of the JFF 
voltage at fixed current on the c-axis magnetic field, Hz 
for the mesa #1. We can see that the JFF voltage is 
highly sensitive to a small component of Hz resulting 
in large negative JFF magnetoresistance Rjff{Hz) = 
U{Hz)/L At the fixed current (0.2 mA) and fixed in- 
plane field (1 tesla) the JFF resistance first decreases 
linearly with Hz and then drops down sharply to the 
low-resistance state. To understand the mechanism of 
this drop, we show in the right panel of Fig. [51 a series of 
I-V dependences traced at the different fixed Hz values 
marked in the left plot. We can see from this plot, that 
the voltage drop to the low-resistance state originates 
from a nonmonotonic character of the I-V dependence 
which appears at high enough Hz. The excess current. 




[Oe] u (V) 



FIG. 5: Left plot shows Hz dependence of the z-axis voltage 
at fixed current / =0.2mA and in-plane filed Hx — IT for the 
mesa #1. The right plot shows the current- voltage (TV) de- 
pendences at several fixed Hz values marked in the left plot. 
The horizontal dashed line marks the current value in the left 
plot. One can see that the I-V dependence acquires a jump 
around Hz ~ 50Oe, which is a coming from the nonmono- 
tonic contribution, SI{U,Hz), from the pancake vortices. As 
a consequence, the Hz (or angular) dependence of voltage at 
fixed current in the left side also has a jump. 



#2Pristine, T = 50K 




-0.8 
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FIG. 6: Plots of the ratio Sj/Hz vs voltage for H^ = IT 
and three tilt angles demonstrating that the excess pancake 
current is proportional to concentration of PVs. 



SI{U,Hz) = I{U,Hz) - /(/7,0), in the I-V dependences 
is a direct consequence of an additional damping of the 
moving JV lattice induced by the PV stacks, as discussed 
in the theoretical sections above. Figure IHl shows plots 
of the ratios Sj/Hz vs bias voltage U for the mesa #2 
before irradiation at different tilt angles corresponding 
to different Hz. One can see nonmonotonic dependence 
5j(U)^ that first increases rapidly, reaches a maximum 
and then decreases slowly. Another remarkable feature 
is proportionality of 5j to Hz up to Hz ~ 87 Oe (0.5° 
tilt), which is seen from an almost ideal collapse of the 
5j /Hz curves at different tilt angles. Figure [7| shows the 
I-V dependences of the mesa #1 for two values of the 
in-plane field, H^ = 1.04T and 2.08T, and for very close 
values of c-axis field, Hz ~ 40.5Oe. We can see that for 
Hx = 2.08T the amplitude of the excess current is more 
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FIG. 7: The current- voltage dependences of the mesa #1 for 
two values of the in-plane field, Hx = 1.04T and 2.08T, and 
for very close values of c-axis field Hz = 40.5 and 40.8Oe. 
For reference, the I-V dependences at Hz = OOe are also 
shown for both Hx- The amplitude of the excess PV current 
for Hx — 2.08T is more than four times smaller than for 
Hx 1.04T. 
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FIG. 8: Temperature evolution of the excess PV current den- 
sity for the mesa #2 at Hx — IT and tilt angle 0.3° cor- 
responding to Hz = 52.40e. The voltage corresponding to 
maximum excess current decreases with decreasing tempera- 
ture. 



than four times smaller than for Hx = 1.04T, which is 
consistent with the theoretical estimates. Figure |H1 illus- 
trates the temperature dependence of the excess current 
for the mesa #2. We can see that at lower temperatures 
maximum in the excess current increases and moves to 
the lower voltages. We also note that for the overdoped 
mesa the excess PV current density is 4-5 times larger 
than for the underdoped one. 

Figure El shows comparison of the PV-induced damp- 
ing for the mesa #2 before and after irradiation for three 
temperatures, 40K, 50K, and 60K. It is interesting to 



FIG. 9: Influence of heavy-ion irradiation on the pancake 
excess current. The cur rent- volt age characteristics are shown 
for mesa #2 before and after irradiation, for the magnetic 
field IT tilted at angle 0° and 0.3° (corresponding to Hz ~ 
520e). Irradiation dose corresponds to the matching field, 
— GOG. Irradiation does not influence motion of the JV 
lattice without PVs, at Hz = but strongly reduces pancake 
effect at T = 40K and 50K. This reduction is much weaker at 
60K. 



note that the columnar defects induced by heavy-ion ir- 
radiation have negligible direct influence on dynamics of 
JV lattice (IVs at 0° tilt angles are not noticeably in- 
fluenced by the irradiation). On the other hand, they 
strongly influence the PV-induced damping at low tem- 
peratures. We can see that at 40K and 50K pinning of 
the PV stacks by the columnar defects almost completely 
eliminates the PV extra damping. At 60 K this damping 
is restored most probably due to the thermal depinning of 
the PV stacks from the columns. This behavior is consis- 
tent with location of the thermal depinning temperature 
T* evaluated in Sec. Ill PI (for underdoped BSCCO we 
estimated T* - 45K). 



IV. DISCUSSION 

Strong dynamic interaction of the JVs with the PV 
stacks leads to a very high sensitivity of the JV trans- 
port to a very small concentration of the PV stacks which 
we observed experimentally. At a given driving force (or 
DC current) that interaction slows down the JV lattice 
and thus reduces the JFF voltage leading to the pro- 
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nounced negative JFF magnetoresistance with respect to 
the c-axis component of the magnetic field. On the other 
hand, at a given voltage the extra damping due to the 
PV stacks results in increase of the current. Theoreti- 
cal description of these effects is based on calculation of 
the oscillating zigzag deformations of the PV stacks in- 
duced by the moving JV lattice and the average friction 
force due to these deformations. Several features of the 
experimental PV excess current, such as its nonmono- 
tonic dependence on voltage, its decrease with increasing 
temperature and the in-plane field, are in a good quali- 
tative agreement with the theory. A strong suppression 
of the PV damping by columnar defects also provides 
confirmation for suggested mechanism. A detailed quan- 
titative comparison between the experiment and theory 
requires an accurate knowledge of the microscopic pa- 
rameters of superconductor, London penetration depth, 
Aa6, anisotropy 7, and pancake fiux-fiow resistivity p^, 
as well as geometrical parameters of the mesas. 

To estimate the key theoretical parameters dp and Vr 
we need values of the magnetic spring constant Kp and 
the PV viscosity coefficient, rjp. The magnetic spring 
constant is mainly determined by the London penetra- 
tion depth. For example, taking a value A = 300nm 
typical for overdoped BSCCO at T = 65K, we obtain 
Kp/s ^ 65N/m^ and Bx ^ 36G. The value of the viscos- 
ity coefficient is more uncertain. To our knowledge no di- 
rect transport measurement of the PV ffux-flow resistiv- 
ity in BSCCO has been published. Microwave measure- 
ments of the ffux-flow resistivity in single crystals (at 40.8 
GHz)2^ and fflms (at 48 GHz)^^ show that r]p ex dp^/dB^ 
has strong temperature dependence and suggest a typical 
value for T = 60K, 7^p/s~10~^-10~^N s/m^ correspond- 
ing to the slope dpf^/dBz = 2-20 pQ. cm/tesla. This gives 
the following estimate for the ffux ffow-conductivity at 
Bz = Bx, crff(5A)-2 • 10^-2 • 10^[n cm]-i. These data 
suggest that the maximum PV effect is expected at very 
small voltage, Vr ~ 0.1-1/iV. Experimentally, the max- 
imum excess current is observed at signiffcantly higher 
voltages, Vr ~ 15-30/iV corresponding to relaxation fre- 
quencies 7-15 GHz. Also, the simple theory does not 
quite describe shapes of the experimental Sj{U) depen- 
dences, dit U > Vr the excess current decreases slower 
than expected 1/U decay. On the other hand, we found 
that shapes of Sj{U) dependences also somewhat differ 
for under doped and overdoped mesas. These discrepan- 
cies can be explained if we (i) assume small- velocity value 
for the viscosity constant, r]p/s ~ 2 • 10~^N s/m^, which 
is signiffcantly larger than the microwave measurements 
suggest, and (ii) assume that it decreases with increas- 
ing PV velocity. In the temperature and in-plane ffeld 
region probed in experiment inffuence of thermal ffuc- 
tuation PV's considered in Sec. IIIBI is signiffcant. For 
example, for the overdoped mesa at 65 K, we estimate 
a J ^ 0.15-0.2 and the reduced temperature given by Eq. 
^asf^ 0.2 for B^ = IT and ^ 0.8 for B^ = 2T. 
Shape deviations, at least partly, can be due to the PV 
thermal ffuctuations. Our calculations show that the PV 



ffuctuations mainly reduce the amplitude of 5j and some- 
what weaken the dependence Sj{U) at large U but they 
do not shift much the maximum of the Sj{U) dependence. 

Another possible reason for disagreement with the sim- 
ple theoretical picture is that at small velocities the 
JV lattice may move inhomogeneously, via thermally- 
activated jumps. This regime of motion is not considered 
by the theory at all and it may lead to an increase of ap- 
parent Vr. We expect that this inhomogeneous regime is 
more pronounced for smaller sample sizes in which av- 
eraging over PV locations is incomplete. Therefore, a 
systematic size dependence of the effect would help to 
understand mechanism of the low- volt age regime. 

Note that, due to a very small value of Vr, without 
scrutinizing LV shape, the low- volt age enhancement of 
damping can be easily misinterpreted as an enhancement 
of the critical current. In our experiments, a detailed 
study of the small- voltage region shows that the JV lat- 
tice slowly moves below the maximum current meaning 
that this maximum is due to the enhancement of damp- 
ing. This, of course, does not exclude possibility that in 
other experimental conditions (e.g., at low temperatures) 
the main PV effect would be enhancement of the critical 
current. 

An extreme sensitivity of the JV ffux-ffow voltage to 
a very small concentration of the PVs suggests a very 
attractive possibility to use small-size mesas for detection 
of penetration of individual PV stacks which will be seen 
as current or voltage jumps. The change of the total 
current due to the penetration event of a single PV stack, 
SIi , is expected to be independent on the mesa size and 
is given by 

- f *»■ 

Its maximum value, ^/i,max, is realized ditV = Vr and at 
low temperatures, where thermal-fluctuation suppression 
is small, can be estimated as 



J J 



ln(A/r^)- 



Taking jj = 500A/cm^, A = 300nm, and h = 2, we 
estimate ^/i,max ^ 0.5/iA. On the other hand, for the 
underdoped mesa #2 at 50K from Fig. we obtain 
(^A,max ^ 0.12/iA consistent with the above theoretical 
estimate. For the overdoped mesa ^/i,max is 4-5 times 
larger. Such currents are detectable experimentally. We 
conclude that this technique can be practically used for 
detection of the individual PV stacks penetration. 
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